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Controlling the Schottky barrier at MoS 2 |metal contacts by inserting a BN monolayer 
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Making a metal contact to the two-dimensional semiconductor M 0 S 2 without creating a Schottky 
barrier is a challenge. Using density functional calculations we show that, although the Schottky 
barrier for electrons obeys the Schottky-Mott rule for high work function (> 4.7 eV) metals, the 
Fermi level is pinned at 0.1-0.3 eV below the conduction band edge of M 0 S 2 for low work function 
metals, due to the metal-MoS 2 interaction. Inserting a boron nitride (BN) monolayer between the 
metal and the M 0 S 2 disrupts this interaction, and restores the M 0 S 2 electronic structure. Moreover, 
a BN layer decreases the metal work function of Co and Ni by ~ 2 eV, and enables a line-up of the 
Fermi level with the M 0 S 2 conduction band. Surface modification by adsorbing a single BN layer is 
a practical method to attain vanishing Schottky barrier heights. 

PACS numbers: 73.30.+y, 73.20.At, 79.60.Jv 


Introduction. Single layers of transition metal dichalco- 
genides (TMDs) such as molybdenite, MoS 2 , can be 
exfoliated through micromechanical cleavage, similar to 
graphene. 1 In contrast to graphene however, a MoS 2 
monolayer is a semiconductor with a sizable band gap 
of 1.8 eV, 2 which has triggered a large interest in 
TMD semiconductor devices. 3-5 Contacting MoS 2 to 
metal electrodes remains a problem, as it tends to pro¬ 
duce unexpectedly high contact barriers and resistances. 
Early photoemission experiments claimed that the Schot¬ 
tky barriers at MoS 2 |metal interfaces obey the ideal 
Schottky-Mott rule, 6 suggesting the possibility to con¬ 
trol the Schottky barrier height (SBH). In particular, the 
SBH for electrons might be reduced to zero using a metal 
with a sufficiently low work function. However, the more 
recent device experiments do not give zero SBHs, neither 
for metals with high work functions, nor for metals with 
low work functions. 7 12 

SBHs of metal contacts with conventional semiconduc¬ 
tors such as Si often only weakly depend on the metal 
species, and the Fermi level is pinned inside the semi¬ 
conductor band gap. 13 Common models used to explain 
Fermi level pinning rely upon having a strong (chemi¬ 
cal) interaction at the metal-semiconductor interface that 
yields a large density of interface states with energies in 
the semiconductor band gap. Unlike Si, MoS 2 has no 
dangling bonds at its surface. Its interaction with metal 
surfaces should therefore be relatively weak, which makes 
it rather unlikely that midgap interface states are formed 
at a high density. Indeed a recent density functional the¬ 
ory (DFT) study claims there is only a weak Fermi level 
pinning at MoS 2 |metal interfaces. 14 That still leaves the 
prospect of a zero SBH using a metal that has a suffi¬ 
ciently low work function. 

In this paper we study the Schottky barriers at 
MoS 2 |metal interfaces by DFT calculations with the ob¬ 
jective of designing a contact with zero SBH. We start 
from a series of metals covering a wide range of work 
functions (3.8-5.8 eV). Unlike previous studies we do not 
so much focus on the chemical interactions with spe¬ 


cific metals, 15-17 but on establishing general rules for 
the SBHs. We show that for clean MoS 2 |metal inter¬ 
faces with high work function metals the Fermi level 
is not pinned in the MoS 2 band gap, and the SBH 
shows Schottky-Mott behavior. 13 This breaks down for 
low work function metals, and the Fermi level gets pinned 
just below the MoS 2 conduction band, leading to a finite 
SBH in the range 0.1-0.3 eV. The metal-MoS 2 interaction 
at the interface perturbs the electronic structure of MoS 2 , 
its conduction bands in particular, creating a density of 
interface states just below the conduction band that pins 
the Fermi level. 

We “unpin” the Fermi level by inserting a h -BN mono- 
layer between the metal surface and MoS 2 . It breaks the 
direct metal-MoS 2 interaction and destroys the interface 
states. Like graphene, MoS 2 is physisorbed on the UBN- 
covered substrate, 18,19 which leaves its electronic struc¬ 
ture nearly unperturbed. Moreover, adsorption of h -BN 
on a metal surface commonly decreases the work function 
considerably. It turns high work function metals such as 
Co and Ni into low work function substrates. The com¬ 
bined effects of breaking the metal-MoS 2 interaction and 
lowering the metal work function yields zero SBHs for 
contacts between /i-BN-covered Co or Ni and MoS 2 . 

Computational details. We use projector augmented 
waves (PAW) as implemented in the VASP code. 20-23 
The MoS 2 |metal interface is modeled by a slab of four 
layers of metal atoms with a monolayer MoS 2 adsorbed 
on one side. We optimize all atomic positions, keeping 
the layer of metal atoms furthest removed from the adsor- 
bant in its bulk geometry. We force the metal lattice to 
be commensurable to the MoS 2 lattice, choosing in-plane 
supercells with a mismatch between the MoS 2 and the 
metal lattices of less than 1%. The metal-adsorbant bind¬ 
ing distance is important for obtaining an accurate inter¬ 
face potential profile, and in some cases this distance de¬ 
pends sensitively on the DFT functional. 24 The PBE gen¬ 
eralized gradient approximation 25 and the optB88-vdW- 
DF van der Waals density functional (vdW-DF) 26,27 give 
MoS 2 |Ag, Au, Pd and Pt(lll) binding distances within 
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0.1 A of one another and interface potential steps within 
0.05 eV. The local density approximation (LDA) 28 yields 
on average ^ 0.2 A shorter binding distances and ~ 0.15 
eV larger potential steps. We use PBE in the following to 
avoid the risk of overbinding commonly found with LDA. 
Whereas the interface potential step does not depend too 
critically on the DFT functional, the work function of a 
clean metal surface Wm can be more sensitive, 14,29 which 
is then reflected in the calculated SBHs, see Eq. (1). 

MoS 2 \metal interfaces. The SBH for electrons can be 
written as 

= AV, (1) 

with Wm the work function of the clean metal surface, 
X the electron affinity of M 0 S 2 , and AV the potential 
step formed at the MoS 2 |metal interface, see Fig. 1. The 
potential step can be calculated without resorting to the 
details of the potential profile across the interface, or its 
electronic structure, because AV = Wm — W a d s |M? where 
Wads| m is the work function of the metal surface covered 
with M 0 S 2 . The results are shown in Fig. 1(a). The 
SBH is then obtained from Eq. (1), using the calculated 
X = 4.30 eV, cf. Fig. 1(b). 

Two regimes can be distinguished in Fig. 1(a). Met¬ 
als with Wm > 4.7 eV yield a similar AV, whereas 
AV strongly depends on Wm for Wm <4.7 eV. The 
cross-over between the two regimes occurs for Cu, where 
Wcu ~ X + AF and the SBH is minimal, see Fig. 1(b). 
The results for the high work function regime are consis¬ 
tent with the common observation that physisorption re¬ 
sults in a net decrease of the work function, i.e., a positive 
AV. This has been explained in terms of the Pauli ex¬ 
change repulsion between the electrons of the metal and 
those of the overlayer yielding a net pushback of electrons 
into the metal. 24 The effect is fairly moderate for M 0 S 2 , 
resulting in AV ~ 0.35 eV for the high work function 
metals, and a size that does not depend critically on the 
details of the metal or the interface structure. 

With a constant AV the SBH simply follows the work 
function, i.e., the slope S = d<fr n /dWM ~ 1, which is 
the Schottky-Mott rule. This rule is typically found in 
the absence of any interface states with energies in the 
semiconductor band gap, which is consistent with the 
MoS 2 |metal interaction being weak. If ideal Schottky- 
Mott behavior would persist for the low work function 
metals, the SBH would vanish for Wm < Wc u - Clearly 
this is not the case in Fig. 1(b). The SBH has a minimum 
at Cu, but it increases again for the low work function 
metals. 

We analyze this behavior for Mg, being the metal with 
the lowest work function in this study. The band struc¬ 
ture of the (1 x 1 ) MoS2|Mg(0001) interface is shown 
in Fig. 2 (a). At the optimized equilibrium distance 
d = 2.2 A, the M0S2 bands are significantly perturbed 
by the interaction with the substrate. This perturba¬ 
tion can be visualized by comparing the density of states 
(DoS) of MoS2|Mg to that of free-standing M0S2, see 
Fig. 2 (e). The MoS2|Mg interaction leads to interface 
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FIG. 1. (Color online) (a) The potential step AV at the 
MoS 2 |metal interface versus Wm, the work function of the 
clean metal surface, (b) The Schottky barrier height (SBH) 
for electrons <f> n . The blue and green dashed lines in (a) and 
(b) indicate the Schottky-Mott rule and Fermi level pinning, 
respectively. Inset: schematic energy diagram of the interface, 
with x the electron affinity of M 0 S 2 . 


states in the M 0 S 2 band gap that are energetically close 
to the bottom of the M 0 S 2 conduction band. This is 
seen most clearly in the difference DoS, i.e. ADoS = 
DoS a ds|M ~ DoS a ds — DoSm, represented by the green 
curve in Fig. 2(e), where DoS a d s |M> DoS a d s , DoSm are 
the DoSs of the interface, the free-standing adsorbate, 
and the clean metal substrate, respectively. 

The ADoS is negligible for energies inside the M 0 S 2 
band gap, except for a region < 0.4 eV below the con¬ 
duction band edge (shaded orange in Fig. 2(e)). The 
negative sign of the calculated potential step at the 
MoS 2 |Mg interface, AV = —0.77 eV, shows that elec¬ 
tronic charge is transferred from the Mg substrate to 
the M 0 S 2 overlayer. These electrons populate the inter¬ 
face states, thereby pinning the Fermi level in the band 
gap. For high work function metals the Fermi level is 
well within the M 0 S 2 gap, where the density of interface 
states is negligible. The corresponding SBHs then obey 
the Schottky-Mott rule. It is only for low work function 
metals, when the Fermi level approaches the bottom of 
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FIG. 2. (Color online) (a) Band structures of the 
MoS2|Mg(0001) interfaces at the equilibrium distance d eq — 
2.2 A and (b) at d = 5 A. The red color measures a projection 
of the wave function on the M 0 S 2 orbitals. The Fermi level 
is set at zero energy. (c,d) The electron difference density 
A n(z), corresponding to the interface distances of (a,b), re¬ 
spectively. (e) The total density of states (red) corresponding 
to (a) DoS ads |M g , (blue) of free-standing M 0 S 2 DoS a ds, and 
(green) the difference ADoS = DoS a d s |M g — DoS a d s — DoSMg. 
The orange shading indicates the position of the gap states, 
created by the interaction at the interface. 


the M 0 S 2 conduction band, that interface states become 
noticeable and pin the Fermi level. 

One can destroy the interface states by breaking the 
M 0 S 2 1Mg interaction. This is demonstrated by Fig. 2(b), 
which gives the band structure of MoS 2 |Mg( 0001 ) with 
the adsorbant placed at an artificially large distance 
d = 5 A from the metal surface. The M 0 S 2 bands are 
unperturbed, as there is no chemical interaction with the 
substrate at this distance. The Fermi level is in the M 0 S 2 
conduction band, which is not surprising as the work 
function of Mg(0001), Wyig = 3.8 eV is much smaller 
than the M 0 S 2 electron affinity x = 4.3 eV. The calcu- 
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FIG. 3. (Color online) Top and side view of the MoS 2 |/i- 
BN|Co(lll) structure. The black rhombus indicates the sur¬ 
face supercell. 


lated interface potential step, AV n [ = —0.5 eV, indeed 
corresponds to Wyig — x, as it should for electron trans¬ 
fer from Mg to M 0 S 2 to equilibrate the Fermi level. It 
results in a zero SBH, according to Eq. (1). 

This interpretation is confirmed by the electron den- 
sity difference A n(z) = n ads | M (2) - n ads {z ) - n M (z), 
where n ads |M(z),^adsM, hm(^) are the plane-averaged 
electron densities of the interface, the free-standing ad¬ 
sorbate, and the clean metal substrate, respectively. At 
an MoS 2 |Mg distance d = 5 A, A n(z) shows an accu¬ 
mulation of electrons at the position of M 0 S 2 , and an 
electron depletion at the Mg surface, see Fig. 2(d), con¬ 
sistent with an electron transfer from Mg to M 0 S 2 , which 
creates an interface dipole and potential step to equili¬ 
brate the Fermi level. At the equilibrium MoS 2 |Mg dis¬ 
tance, A n(z) shows a much more complicated pattern, 
see Fig. 2(c), which is consistent with an interface inter¬ 
action that alters the electronic structure. 

M 0 S 2 1 h-BN\ metal interfaces. A Schottky barrier at a 
MoS 2 |metal contact is unavoidable, as the interface inter¬ 
action leads to states that pin the Fermi level below the 
conduction band. Breaking this interaction by introduc¬ 
ing a vacuum spacing at the interface, is not a practical 
way to make a metal-semiconductor contact. However, 
inserting an inert layer between the metal and M 0 S 2 can 
be. This layer has to be sufficiently thin as not to form 
a high barrier for electron transport. A monolayer of 
h -BN is ideally suited. A single layer of h -BN can be 
deposited or grown on a range of metal substrates, and 
it is stable under ambient conditions. Substrates con¬ 
sisting of a transition metal (111) surface covered by a 
h -BN monolayer, are widely available. Sandwiching a 
monolayer of h -BN between two metal electrodes gives 
metallic conduction, 30 which indicates that the layer is 
transparent to electrons. 

We build MoS 2 |/i-BN| metal structures by putting a 
(2 x 2) M 0 S 2 cell on top of a (y/7 x y/7) /i-BN|Co or 
Ni cell, see Fig. 3, with the h -BN layer in a (1 x 1) reg¬ 
istry with the underlying metal surface, as in Refs. 24 
and 29. As before, we fix the lattice constant of M 0 S 2 , 
and adapt the lattices of h -BN and the metal(lll) sub- 
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FIG. 4. (Color online) AV versus Wm■ Inset: 4> n versus Wm- 
The blue and green dashed lines indicate the Schottky-Mott 
rule and Fermi level pinning, as in Fig. 1. 


strates accordingly, which requires a 4% squeeze of the 
h -BN lattice. A h -BN monolayer is chemisorbed on Co 
and Ni(lll) surfaces, but M 0 S 2 and h -BN are bonded 
by a weak, van der Waals, interaction. Such interactions 
are not represented in the PBE functional, so we use the 
optB88-vdW-DF van der Waals density functional here. 

Adsorption of h -BN has a dramatic effect on the work 
function; it reduces the work functions of Co and Ni(lll) 
by 1.9 eV and 1.8 eV, respectively. These reductions re¬ 
sult from large interface dipoles that are formed at the 
/z-BN| metal interfaces, where Pauli exchange repulsion 
between the electrons at the interface gives an important 
contribution. 24 Adsorbants in the form of self-assembled 
monolayers (SAMs) are commonly proposed in order to 
modify substrate work functions. However, SAMs often 
suffer from disorder, which diminishes their effect. Ad¬ 
sorption of h -BN leads to a well-defined structure that 
is much less susceptible to disorder, and gives a sizable 
work function lowering. Moreover h -BN presents a sur¬ 
face that not only is chemically relatively inert, but also 
does not change its structure upon adsorbing further lay¬ 
ers. 

The implications of work function lowering by h -BN 
adsorption are clearly demonstrated in Fig. 4. Direct 
adsorption of M 0 S 2 on Co and Ni(lll) gives a behavior 
that is typical for high work function metals. The poten¬ 
tial step at the MoS 2 |metal interface is AV ~ 0.35 eV, 
and the SBHs follow the Schottky-Mott rule. In contrast, 
adsorption of M 0 S 2 on ft-BN|Co and Ni(lll) substrates 
gives a negative AF, and it gives a zero SBH. Insert¬ 
ing a h -BN layer has not only effectively decreased the 
substrate work function, but it has also weakened the 
M 0 S 2 1 substrate interaction that yielded Fermi level pin¬ 
ning and nonzero SBHs for clean low work function metal 
substrates. 




FIG. 5. (Color online) (a) Band structures of the 
MoS2|Co( 111) and (b) the MoS 2 |/i-BN|Co(lll) interfaces. 
The red color measures a projection of the wave function on 
the M 0 S 2 orbitals. The Fermi level is set at zero energy. (c,d) 
The corresponding electron difference density A n(z). The or¬ 
ange line indicates the position of the /z-BN layer. 


The difference between adsorbing M 0 S 2 directly onto 
a clean Co surface and onto a h -BN covered Co sur¬ 
face is observed in the corresponding electronic struc¬ 
tures shown in Fig. 5. Direct adsorption onto Co(lll) 
perturbs the M 0 S 2 bands considerably, due to the inter¬ 
action at the interface, as shown in Fig. 5(a). In contrast, 
adsorbing M 0 S 2 onto a /z-BN|Co(lll) substrate hardly 
perturbs the M 0 S 2 bands at all, as demonstrated by Fig. 
5(b), which is a clear indication that the interaction be¬ 
tween h -BN and M 0 S 2 is weak. 

Adsorbing M 0 S 2 directly onto Co leads to the Fermi 
level being well within the M 0 S 2 band gap, which is not 
surprising as Co is a high work function metal. In con¬ 
trast, adsorbing M 0 S 2 onto a ft-BN|Co substrate yields 
a Fermi level that is pinned at the bottom of the M 0 S 2 
conduction band. The electron density differences A n(z) 
shown in Fig. 5(c) and 5(d) substantiate this picture. 
Figure 5(c) shows the complicated pattern that is typ¬ 
ically associated with a direct M 0 S 2 Imetal interaction, 
compare to Fig. 2(c). A n(z) for M 0 S 2 on a /z-BN|Co 
shows an accumulation of electrons at the position of 
M 0 S 2 , and an electron depletion at the position of the 
h- BN, which reflects an electron transfer from the sub¬ 
strate to M 0 S 2 to equilibrate the Fermi level, compare 
Figs 5(d) and 2(d). 

Conclusions. We have shown that contacting M0S2 
with low work function metals leads to Fermi level pin- 
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ning at 0.1-0.3 eV below the conduction band edge. This 
behavior results from the interaction at the interface be¬ 
tween the metal and the M 0 S 2 , creating a considerable 
density of interface states just below the M 0 S 2 conduc¬ 
tion band. Inserting a boron nitride (BN) monolayer 
between the metal and the M 0 S 2 destroys these interface 
states, and recovers the unperturbed M 0 S 2 band struc¬ 
ture. In addition, absorbing h -BN on Co(lll) or Ni(lll) 
decreases the metal work function by close to 2 eV. We 


predict that contacting M 0 S 2 with ft-BN|Co or Ni(lll) 
does not give a Schottky barrier. 

Acknowledgement. We thank Menno Bokdam for use¬ 
ful discussions. This work is part of the research program 
of the Foundation for Fundamental Research on Matter 
(FOM), which is part of the Netherlands Organisation for 
Scientific Research (NWO). The use of supercomputer fa¬ 
cilities was sponsored by the Physical Sciences Division 
(EW) of NWO. 


* g.h.l.a.brocks@utwente.nl 

1 A. K. Geim and I. V. Grigorieva, Nature 499, 419 (2013). 

2 K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, 
Phys. Rev. Lett. 105, 136805 (2010). 

3 B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 
and A. Kis, Nat. Nanotechnol. 6, 147 (2011). 

4 Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, 
and M. S. Strano, Nat. Nanotechnol. 7, 699 (2012). 

5 H. Wang, L. Yu, Y.-H. Lee, Y. Shi, A. Hsu, M. L. Chin, 

L. -J. Li, M. Dubey, J. Kong, and T. Palacios, Nano Lett. 
12, 4674 (2012). 

6 J. R. Lince, D. J. Carre, and P. D. Fleischauer, Phys. Rev. 
B 36, 1647 (1987). 

7 H. Liu, A. T. Neal, and P. D. Ye, ACS Nano 10, 8563 

( 2012 ). 

8 S. Das, H.-Y. Chen, A. V. Penumatcha, and J. Appen- 
zeller, Nano Lett. 13, 100 (2013). 

9 J.-R. Chen, P. M. Odenthal, A. G. Swartz, G. C. Floyd, 
K. Y. L. Hua Wen, and R. K. Kawakami, Nano Lett. 13, 
3106 (2013). 

10 M. Fontana, T. Deppe, A. K. Boyd, M. Rinzan, A. Y. Liu, 

M. Paranjape, and P. Barbara, Sci. Rep. 3, 1634 (2013). 

11 N. Kaushik, A. Nipane, F. Basheer, S. Dubey, S. Grover, 
M. M. Deshmukh, and S. Lodha, Appl. Phys. Lett. 105, 
113505 (2014). 

12 J. Kang, W. Liu, and K. Banerjee, Appl. Phys. Lett. 104, 
093106 (2014). 

13 R. T. Tung, Appl. Phys. Rep. 1, 011304 (2014). 

14 C. Gong, L. Colombo, R. M. Wallace, and K. Cho, Nano 


Lett. 14, 1714 (2014). 

15 I. Popov, G. Seifert, and D. Tomanek, Phys. Rev. Lett. 
108, 156802 (2012). 

16 W. Chen, E. J. G. Santos, W. Zhu, E. Kaxiras, and 
Z. Zhang, Nano Lett. 13, 509 (2013). 

17 J. Kang, W. Liu, D. Sarkar, D. Jena, and K. Banerjee, 
Phys. Rev. X 4, 031005 (2014). 

18 M. Bokdam, P. A. Khomyakov, G. Brocks, Z. Zhong, and 
P. J. Kelly, Nano Lett. 11, 4631 (2011). 

19 R. Gillen, J. Robertson, and J. Maultzsch, Phys. Rev. B 
90, 075437 (2014). 

20 G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 

21 P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

22 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

23 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

24 M. Bokdam, G. Brocks, and P. J. Kelly, Phys. Rev. B 90, 
201411(R) (2014). 

25 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

26 M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and 
B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004). 

27 J. Klimes, D. R. Bowler, and A. Michaelides, Phys. Rev. 
B 83, 195131 (2011). 

28 J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

29 M. Bokdam, G. Brocks, M. I. Katsnelson, and P. J. Kelly, 
Phys. Rev. B 90, 085415 (2014). 

30 O. Yazyev and A. Pasquarello, Phys. Rev. B 80, 035408 
(2009). 



